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Abstract 

The Casimir energy for a massless scalar field between the closely 
spaced two concentric D— dimensional (for D > 3) spheres is calculated by 
using the mode summation with contour integration in the complex plane 
of eigenfrequencies and the generalized Abel-Plana formula for evenly 
spaced eigenfrequency at large argument. The sign of the Casimir en- 
ergy between closely spaced two concentric D-dimensional spheres for a 
massless scalar field satisfying the Dirichlet boundary conditions is strictly 
negative. The Casimir energy between D — 1 dimensional surfaces close 
to each other is regarded as interesting both by itself and as the key to 
describing of stability of the attractive Casimir force. 

PACS number(s): 03.70.+k, ll.10.Kk, ll.10.Gh, 03.65.Ge 

1 INTRODUCTION 

The measurable consequences of the macroscopic phenomena in the quan- 
tum theory of fields is Casimir effect. This effect due to the vacuum polariza- 
tion of the quantized field was originally derived by H. B. G. Casimir pQ. He 
calculated the negative renormalized quantum vacuum energy for the electro- 
magnetic field bounded by two uncharged parallel plates. He concluded that 
there must exist an attractive force between the plates. The attractive force: 
F = — ^Qai which is expressed in terms of the Planck constant H, the velocity 
of light c, and the distance between the plates a, should act on unit area of two 
uncharged conducting plane parallel plates in vacuum. This attractive force was 
confirmed experimentally [2l [3J |4] . Both theoretical and experimental studies 
of the Casimir effect can provide great insight into understanding the nature of 
the quantum vacuum |H HI 13 HI IS] 

At first time, having found the negative energy for uncharged conducting 
parallel plates due to Casimir effect, the hope was what the same would appear 
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for the spherical geometry. This expectation was shattered by Boyer [TU] . Boyer 
first showed Casimir energy for the spherical shell is positive. The sign of the 
positive Casimir energy produces the repulsive force. Boyer's result has been 
later confirmed by using different regularization techniques [TTJ [T2J [T3J . Nowa- 
days, the nature of the Casimir energy is the strong dependence on the geometry 
of the spacetime, the dimension of the spacetime and on the boundary condition 
imposed [7l ll4[[T5] . Recently, the Casimir energy between two concentric spheres 
and cylinders in D = 3 dimension have been considered by using the different 
regularization methods [ig[]7|[IEl[ISI13l2Ill21ll2l21l^ 
for the scalar and electromagnetic fields. Moreover, the investigation of the di- 
mensional dependence of the Casimir energy is of interest. The Casimir energy 
in the higher dimensional spacetime has a long history. One of the pioneering 
works is that the Casimir energy for a massless scalar field and a massless vec- 
tor field in a D— dimensional rectangular cavity were derived by Ambjorn and 
Wolfram [35]. Afterwards, the scalar and electromagnetic Casimir energy for 
a D— dimensional sphere were calculated using Green's functions method and 
zeta function regularization [33J 132 133 [351 [37] , and the Casimir interaction be- 
tween two concentric spheres in Z?-dimensions for the scalar fields with Robin 
boundary condition and the electromagnetic fields at finite temperature have 
been studied [35J [3J5] . In particular today, physical theories and models with 
extra dimensions are active areas of research, for example String theory and 
Brane world models. 

As matter of fact, the sign and magnitude of the Casimir energy may 
strongly depend on (a) the spacetime dimensionality, (b) the type of the bound- 
ary conditions, (c) type of the fields, (d) the lengths between the surfaces is crit- 
ically dependent on their nanometre-scale shape, (e) the curved spacetime back- 
ground, (f) compactness spacetime, and (g) the finite temperature [7JHU]. In 
this work, the consequences of (c), (e) and (g) will not be considered, and we con- 
sider the Casimir energy between closely spaced two concentric D— dimensional 
spheres for a massless scalar field satisfying the Dirichlet boundary conditions. 
From a mathematical point of view, a simple method for calculating the Casimir 
energy between closely spaced two concentric D— dimensional spheres for a mass- 
less scalar field is developed which is based on a direct mode summation with 
the contour integration in a complex plane of eigenfrequencies and using the 
generalized Abel-Plana sum formula for evenly spaced eigenfrequency at large 
argument. One of the motivations of our work for us to perform the Casimir 
energy calculation for a massless scalar field between D — 1 dimensional surfaces 
close to each other is that the physics in higher dimensional spacetime have be- 
come a trend since the existence of the extra dimensions allows the solving of 
the some fundamental problems in physics as the hirarchy problem. 

The organization of this paper is as follows. In section 2, The Casimir 
energy of a massless scalar field subjected to the Dirichlet boundary condi- 
tions on between closely spaced two concentric D— dimensional spheres is cal- 
culated without any approximation techniques. Concluding remarks and dis- 
cussion of the Casimir energy for a massless scalar field in an annular region of 
D— dimensional geometry is presented in section 3. 
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The units are such that h = c= 1 . 



2 CASIMIR ENERGY 

We begin with a massless scalar field considered in a D— dimensional spher- 
ical geometry, where the metric is given by 

ds 2 = dt 2 — da 2 , where 

da 2 = dx\ + dx\ + dxj + + dx 2 D . (1) 

In spherical coordinates 

x\ — r sin#i sin02 sin03 sin0£>_2 cos</) 

x 2 = r sin#i sin0 2 sin0 3 sin0£>^ 2 sin</> 



XD-l 

x D 



r sin cos 6*2 
r cos 0\ . 



(2) 



Where < 6j < vr, j = 1, 2, 3, D-2 and < < 2tt. a;?+a:|+a;§+. 
r 2 is defined by D dimensional hypersphere. The metric becomes 



ds 2 = dt 2 -[dr 2 + r 2 dQ\ + r 2 sin 2 6»id6»| + .... + r 2 sin 2 1 sin 2 2 sin 2 D - 2 d(j) 2 ] . 

(3) 

Massless scalar field satisfies the Klein Gordon equation in this geometry is 
given by 

□* (t, r, 0, 4>) = (where = 0i,02, ,^-2). (4) 

Here □ is the D'Alembertian operator associated with the metric given by the 
line element in Eq. (3). Solution of Eq. (4) could be easily found by using the 
method of separation of variables and is written as (for D > 4) 



tf(t,r,t?,0) = ^ [Akr-^JA^+Bur-^Nviur) 

{A} 

j [] sin M — 0, cZt::tS^ 2 (cos0,)| (0 C _ 2 ,, 



(5) 
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Where v = k + and {A} refers to m, I (= M x ) , M 2 , M 3 , M 4 , , M D _ 3 , 

and (= Md_ 2 ) . D is the number of dimension of sphere. J v (wr) and N v (cor) 
are Bessel functions first and second kind, respectively. And the function 

^m D -"-i-m d 1 _ 2 ( cos ^m) corresponds to the Gegenbauer or ultraspherical poly- 
nomials. Yf m (6d-2, (f>) are the spherical harmonics. We note that the spherical 
boundary condition at & (t,r = a, "&,4>) = and (t, r = 6, 0) =0 (a (6)is 
inner (outer) radius of D-dimensional sphere) has not imposed on Eq. (5) yet. 
Hence, uj still remains a continous parameter, while m, f (= Mi) , M 2 , M 3 , k (= Md_ 
take the values 

m = ,0,1,2, , 

^(=Mi) = 0,1,2,3,4, ,M 2 

M 2 = 0,1,2,3,4, ,M 3 

M 3 = 0,1,2,3,4, ,M 4 



Md-4 = 0,1,2,3,4, ,M D - 3 

M D _ 3 = 0,1,2,3,4, ,M D _ 2 

k(=M D _ 2 ) = 0,1,2,3,4, (6) 

We now impose the boundary conditions for D— dimensional spherical ge- 
ometry i.e., 

^uim (t, r = a,'d,(t>) = 0, and similarly (t, r = b, i?, </>) = . (7) 

The eigenfunction that satisfy the boundary conditions is 

(P-2) 



{A} 



(8) 



Where c is the normalization constant and uj is the root of the following tran- 
scendental equation 

J v {ujb) N u (ujo) - J v {uja) N u (ujb) = where v = k + ^—7-—^. (9) 



We define the Casimir energy between two concentric Z?-dimensional spheres 
for a massless scalar field 
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Ec = 2 E WA ' 

{A} 

oo (d-3 Ai/^ + i "I £ oo 

= E nE EE — 

1 oo oo 

= E^M E^ 

(10) 

Where w„£ are eigenfrequencies stems from root of the transcendental equation 

given in Eq. (9) and the degenaracy of each eigenfrequency g( D > (y) (D represents space dimension) 

could be written as 



9^ M = 2v, [V :? {D -^ „ and „ = k + ^ . (11) 
y v ; (D-2)![i/-|(D-2)]! 2 v y 

The zeros of the frequency equation are real and simple since v is real and 
a and b positive. We know that Bessel's series equation are convergent for all 
values of the argument. When |w| is very large the convergence is so slow. 
To render the series useless for the frequency calculation we need the rapidly 
convergent evaluation of the Bessel's function formula. A very rapidly conver- 
gent evaluation of the frequency equation can be obtained by using the uniform 
asymptotic expansions for the Bessel function. Hence, we should examine the 
behavior of the eigenfrequency spectrum for large arguments at fixed k and large 
order as k — > oo. Thus, to carry out the summation with respect to k in Eq, 
the sum Y^Li u nk given in equation (10) replaced by Y^=i ^nfc + X^Li <^nk 
where ui n k is the eigenvalue spectrum of the limit w — > oo at fixed k [101 UHl US] . 
Then, the Casimir energy which is defined by the eigenfrequency spectrum for 
large arguments at fixed k and large order as k — > oo and can be written as 



E C = \ E 9 iD) M E ^k + E «nk ■ (12) 
k=0 \n=l n=l ) 

Now, we calculate the eigenfrequencies for large arguments at fixed v. We em- 
ployed the Hankel's asymptotic expansion [TH HH [41] when v is fixed, wa> 1 
and wt> 1, one obtains 



J v (wa) 



/ 2 /„ V 7T\ Uv ~ l) . /„ V 7T\ 

\ / — cos wa 7r — — sm ua n 



(13) 
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cos [ua - -* - - ) + sin [ua - - - ) 

(14) 

and similar expressions for J v (cob) and N v (ujb) with a interchanged for b. 
Putting (13) and (14) in the frequency equation given by (9), we obtain the 
zeros of frequency equation are almost evenly spaced for very large argument at 
fixed v. 

— ] +% where re = 1,2,3,4,5,6, ... . (15) 

— a J ab 

And, the frequency equation ZJ n k the first sum is given in Eq. (12) as the 
uniform asymptotic expansion of the Bessel function can be written as 

fvlvui, A) = J u (vUj) N v (iscuX) — J v (iaUA) N v (fuJ) (16) 

Where A = f (b > a) and v = k + { -^-. 

Then, the Casimir energy between the closely spaced two concentric D- 
dimensional spheres for a massless scalar field which is defined by the eigenfre- 
quency spectrum for large arguments at fixed k and large order as k — > oo can 
be written as [HI US] 



2 

N v (wo) 




1 OO OO 




2 ^ 

k=0 

= E { c D) +E^. (17) 

Where ZJ n k is the root of the frequency equation given in Eq. (16). 

We consider the first sum defined in Eq. (17). This divergent expression can 
be rendered finite by the use of a cutoff or convergence factor. Then we define 
the first sum, 



e ( c d) - JEs (0) m E^^ 1 

fc=0 71=1 
1 OO 

= \J2 9 {D) M S k , (18) 

where u n k refers to ZJ n k and the factor of e~ auirlk plays the role of an exponential 
regulator which effectively suppresses the high frequency contributions to the 
Casimir energy, and Sk — Y^ n =l Wnfc e~ ctWnk is generated by the frequency 
equation (16). To evaluate the sum Sk, we use the integral representation 
from the Cauchy's theorem |T3j [TBI US [26] that for two functions fk{z) and 
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4>{z) analytic within a closed contour C in which fk(z) has isolated zeros at 

&1,&2 i ^3, i 

^£«M(*)^W*(*)=I>(*i)- (19) 
We choose <f>{z) = z e~ az where a is a real positive constant thus leads to 

±£ dz e-- Z ±ln fe ( Z )=Y /Zj e-^. (20) 

Using this result to replace the sum Sk by a contour integral, the first term of 
the Casimir energy becomes 

Ec = \jt gW(u) ±<f dze- a *z-fhif v (vz,\), (21) 
fc=0 m Jc z 

where the frequency function f v (vz,X) is given Eq. (16). The contour C en- 
closes all the positive roots of the equation f v (uz,X) = 0. This contour can 
be conveniently broken into three parts [13j [TBI [19] 26 . These consist of a cir- 
cular segment Cr and two straight line segments Ti and forming an angle 
4> and 7r — <f> with respect to the imaginary axis. When the radius V is fixed, 
the contour Cr encloses a finite number of roots of the equation f v {vz, A) = 0. 
Since the sum of these roots is obviously infinite, the radius T is a regularization 
parameter, and taking the limit T — > 00 ( when a > ) means the removal of 
the regularization, the contribution of Cr vanishes provided that 0^0. Hence 
the exponential regulator in the Cauchy integral plays the role of the eliminate 
of the contribution to the circular part of the contour integral. Taking the 
contributions along Ti and T2 which are complex conjugates of each other and 
rescaling of integration variable, then Eq. (18) becomes [T2l IT5] 



1 r°° r\ 

E [ c } = - J- hm V (u) Re e"* / dy e - ia ^/b ln f v ( uye -i^ A ) . 
2tt6 a-H) J Q dy 

(22) 

Now we use the Lommel's expansions or the multiplication theorem for the func- 
tion of /j/(ye _i< ^, A) and uniform asymptotic expansion of the modified Bcsscl 
functions one obtains [T51 Ql|l HI] 
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y -In f v {yye-^,X) 



..a (1-A 2 ) 3 
1 12 ^ (^e-^) +^ — 

(1-A 2 ) 4 



24 



+ 



+ 



+ 





720 


(1 


-A 2 ) 5 




1440 


(1 


-A 2 ) 6 



{vye-^f (-v 2 + 19) - {vye~^) A 



-3 {vye^f (v 2 - 9) - 2 (z/t/e"^) 4 



120960 



(i^ye-^) 2 (4i/ 4 - 290z/ 2 + 1726) + (i/ye - **) 4 (8^ 2 - 149) 
+ 4 (vye'^) + [Terms in even powers of (yye , (23) 



where |A 2 - l| < 1 and v = k + ^j^. Inserting Eq. (23) into Eq. (22) 
and using the following integral result 



I(2n) = e 



■f 

Jo 



' b \ 2 " +1 

• (_)"+! ( 2n )! ( -J , where n = 0,1, 2, 3... (24) 



then Eq. (22) becomes 



1 



27r6 q^o 



lim J2 9 (D) { v ) Re 



(i-A 2 y 



fe=0 



12 U 



+i 



+i 



(1 


-A 2 ) 4 




720 


(1 


-A 2 ) 5 




720 


(1 


-A 2 ) 6 



2^ 2 (-^ 2 + 19) (-) + 24^ (- 
a I \a 



■*V-9)(i) + 48,<(A 



+i ^ '— ^(8j/- 580^ 2 + 3452) - - 24is 4 (8v 2 - 149) 

120960 v ; V"/ 

[Terms in imaginary number and even powers of v]} . 



(- 



+ 2880i/> - 



b_V] 

(25) 



All terms in the above equation have the singular term in the regulator param- 
eter a and purely imaginary. Taking the real part of the parenthesis, thus it 
leaves the zero result. 
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(26) 



The meaning of this result is that there is no contribution from a large order as 
k — > oo modes for the Casimir energy between the closely spaced two concentric 
D— dimensional spheres. 

Thus, D space dimension of the Casimir energy given in Eq. (17) included 
high cigcnfrcquency modes i.e. uj — > oo at fixed k can be written as 

fc=0 n=l » v 7 

Where f = + ^r 2 ^ . This divergent sum can be regularized by using the 
Abel-Plana sum formula which could be given as [5J |5] 



Reg 



oo 



n=0 



/(«) - /(-it) 



i 



dt. 



(28) 



Where / (z) is an analytic function in the right half plane and Reg refers to the 
regularized value of the sum. The other useful Abel-Plana sum formula for the 
half integer number is 



/(«) - /(-it) 



i 



dt. 



(29) 



We can rewrite the sum given in Eq. (27) using the Abel-Plana sum formula, 
which leads to 



-. oo 

E [ c D) = W 9 l °H») ^g 



fc=0 



E 



b — a 



ah 



j= h E - m - ^ E - 2 M / [y 2 - 1] 1 • 

v/afe fe=Q 2nVab fe=Q 7 c* v - 1 

(30) 



fe=0 



where £ = and d — b ~ a. The first divergent sum in Eq. (30) can be 
removed by using the Hurwitz zeta function 43 
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E * 



fe=0 



P-2)V /_ QD-2) 

H P ' 2 



p + 1 



(31) 



Where p = 0, 1, 2, 3, .. and B p+ i ( — ' 2 2 - ) j is the Bernoulli polynomial. To regu- 
larized the second term in Eq. (30) we use the half integer Abel-Plana sum for- 
mula given in Eq. (29) for D— odd space dimension and the Abel-Plana sum for- 
mula given in Eq. (28) for D— even space dimension, thus the Casimir energy per 
unit surface area on the inner sphere for the specialized cases where our space has 
dimensions D — 4,5, 6, 7, 8, 9, 10, 11 ( for D = 3 result is given by in [19] ) can be 

written as (the total spherical surface area of D dimension is A™ =2^^a^ 



A w 



C _ 

(5) 



1287T 2 



32tt 3 



ab\ £(5) 
d 4 ' 



ah \ C (6) 
If 



1 2 C(4) 



8' C(6) 

l^CCf) _ i 4 C(2)- 

64 ' C(6) 96 ' C(6). 



A" 



(6) 



15 
1024tt 3 



ab\ C(7) 
<i 6 



1-^ 



:C(5) 



15' C(7) 



10 



E. 



(7) 



E. 



(8) 

c 



(9) 



E, 



c 



A" 



ab\ C(8) 



1287T 4 



155 



'15367T 77 C(8) ^ 768tt 2 '' C(8) ^ 256'' C (8) 

7 



>C(6) 



d 7 
35 



1 77 



,C(6) 



12'' C(8) " r 64 7/ C(8) 



; C(4) 



105 

81927T 4 

— ( 

1287T 5 



64 4 C(5) 



7 ' C(9) 525 ' C(9) 



175 , C(4) 



ab\ C (10) 
4445 



!_35 aC(8) 



259 4 C(6) 
77 



7168^ C (10) 491527T 6 '' £(10) 73729'' C (10) 



48 ' C(10) 1152 ' £(10) 
C(8) 5425 8 C(6) 



27195 



C(4) 



175 ^£(2) 



11059207T 2 ' C (10) 86016' C (10) 



E, 



(10) 



.4 1 



945 



65536tt 5 



ab\ C (11) 
d 10 



256 6 C(5) 



3675~ ?? C (11) 



C(9) , 16 4 C(7) 



C(ll) 45' £(11) 



E, 



(ii) 



A 



<ii) 



15 



5127T 6 

3229 



10 

ah \ C(12) 
d 11 

245 



C(6) 



■ 21 2 £(10) 329 
8 C(4) 



23040 ' C (12) 16384 ' £(12) 



10731 
1310727T 8 

10199 
196608tt 4 

22603 
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cog) 

C(12) 

io C(6) 
C(12) 

io C(4) 



6223 



,io 



C(8) 



65536tt 6 ' C (12) 



569625 



14745607T 2 "' 7 C (12) ' 8388608'' C (12) 



C(2) 



4 C(8) 



640 ' C (12) 



(32) 



Where £ (s) is the Riemann zeta function and ri — -p=. Our results are 

V ab 

interest at the limiting case which is narrow slit is defined by r/ — -C 1 [16] . 
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We easily analysis that in the limit a, b — > oo and rf — (77 — ^ and — > 1) 
which means that the surfaces between two spheres converted to the parallel 
plate geometry. One finds that the leading term of the Casimir energy per unit 
area for D dimension can be written as 



E ( C D) 1 r f(D + l) 

A(D) (4.)^ iP 



C(D + 1). (33) 



This result is exactly the same as the scalar Casimir energy of the parallel 
plates for D dimension [32 . Thus our approach developed here has been the 
satisfactory check. As far as we know this result is obtained here for the first 
time. 



3 CONCLUSION 



In this work, we have calculated the Casimir energy between closely spaced 
two concentric D— dimensional spheres for a massless scalar field satisfying the 
Dirichlet boundary conditions. We obtain the numerical results of the Casimir 
energy between the closely spaced two concentric spheres in space dimension 
D = 4 up to D — 11. Although the sign and magnitude of the Casimir energy 
for the spherical shell and cavity change dramatically with the dimension, the 
sign of the Casimir energy between closely spaced two concentric D— dimensional 
spheres does not change. We observed that all spacctime dimensions give us the 
negative renormalized vacuum energy by quantum fluctuations between closely 
spaced two concentric D— dimensional spheres. This result produces to the sign 
of stabilazition of the Casimir energy. The condition of stability will be satisfied 
between D — 1 dimensional surfaces close to each other. 

The interesting result of our calculations is that any approximation tech- 
nique is not needed for our geometry. All contributions in the Casimir energy 
for a massless scalar field comes from the higher frequencies for fixed k between 
two surfaces boundary conditions, k — > 00 frequency modes contribution in 
the Casimir energy is zero for arbitrary width of annular region between closely 
spaced surfaces. Moreover, we find that the Casimir energy per unit area for a 
massless scalar field satisfying the Dirichlet boundary conditions between closely 
spaced two concentric D— dimensional spheres is the same D— dimensional par- 
allel plates in the limit case [32] . 
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